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ABSTRACT 



^ ! We investigate the formation of planetesimals via the gravitational instabihty 

if^ , of sohds that have settled to the midplane of a circumstellar disk. Vertical shear 

^ ! between the gas and a subdisk of solids induces turbulent mixing which inhibits 

'. gravitational instability. Working in the limit of small, well-coupled particles, we 



find that the mixing becomes ineffective when the surface density ratio of solids 
to gas exceeds a critical value. Solids in excess of this precipitation limit can un- 
dergo midplane gravitational instability and form planetesimals. However, this 
saturation effect typically requires increasing the local ratio of solid to gaseous 
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! surface density by factors of two to ten times cosmic abundances, depending on 

the exact properties of the gas disk. We discuss existing astrophysical mecha- 
^ I nisms for augmenting the ratio of solids to gas in protoplanetary disks by such 

^ [ factors, and investigate a particular process that depends on the radial variations 

of orbital drift speeds induced by gas drag. This mechanism can concentrate mil- 
limeter sized chondrules to the supercritical surface density in < few x 10^ years, 
a suggestive timescale for the disappearance of dusty disks in T Tauri stars. We 
discuss the relevance of our results to some outstanding puzzles in planet forma- 
tion theory: the size of the observed solar system, and the rapid type I migration 
of Earth mass bodies. 

Subject headings: planetary systems: formation — planetary systems: proto- 
planetary disks — instabilities — turbulence 
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1. Introduction 

The planetesimal hypothesis states that terrestrial planets and the icy cores of gas giants 
formed in a disk by the accretion of smaller solid bodies, called planetesimals. Planetesimals 
are defined as primitive solids of kilometer size or larger. Such a scale is meaningful because 
it delimits where pairwise gravity dominates gas drag, leading to the runaway or oligarchic 
growth of protoplanets (Lissauer 1993). Because planetesimal formation must occur, at least 
for the giant planets, while gas is still present in the solar nebula, the coupling of gas and 
solids may be crucial to understanding the process. 

The prevailing view for planetesimal formation posits agglomerative growth from sub- 
micron-sized "interstellar" grains to km-sized bodies (cf. the review of Lissauer 1993). How- 
ever, building planetesimals through pure solid-state sticking forces has many problems. 

First, in the inner solar system, interior to the "snow line," the only solids available 
for the formation of the terrestrial planets were rocks. Everyday experience tells us that 
dry silicate particulates of millimeter and larger size, such as sand, do not stick at almost 
any speed of attempted assemblage. Measurements of interparticle collisions in microgravity 
experiments reinforce these intuitive impressions; indeed, /im sized dust particles disrupt 
larger aggregates upon collision at relative velocities > m/s (Blum & Wurm 2000). For 
collisions between fluffy mm-sized agglomerations of fim sized grains, disruption also occured 
at speeds > m/s. However at lower speeds, down to 15 cm/s, only restitution (bouncing) 
was observed, without any sticking (Blum & Muench 1993). Relative velocities considerably 
higher than 1 m/s arise from chaotic motions in a turbulent disk or from differential orbital 
drift in a laminar disk (WeidenschiUing & Cuzzi 1993). For mm-sized crystaUine silicates, 
such as chondrules (see below) , significant shattering occurs only at collision speeds > several 
km/s (Jones, Tielens, & Hollenbach 1996). Thus, compact pieces of rock that make up the 
bulk of the material of chondritic meteorites, whose parent bodies are primitive asteroids 
(rocky planetesimals), are unlikely either to agglomerate or to fragment at the collisional 
velocities common in the nebular disk when there is still appreciable gas present to exert 
appreciable drag on the solids (the well-coupled limit, see Appendix). 

Second, at temperatures significantly below their melting points, even ices in the outer 
solar system may not be much more sticky. Experiments by Supulver et al. (1997) found 
that water frost has spring-like properties and can only induce sticking for collision speeds 
< 0.5 cm/s. Furthermore, if ice balls could grow by continued agglomeration until disruptive 
tidal forces became stronger than the cohesive strength of crystaUine ice, we should expect 
many icy particulates in Saturn's rings to acquire sizes of order ~ 10 — 100 km. In fact, 
except for the occasional embedded moonlet (whose origin may lie in the fragmentation of 
yet larger bodies rather than from the assemblage of ordinary ring particles) , the particulates 
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in Saturn's rings have a maximum size ~ 5 m (Zebker, Marouf, & Tyler 1985), intriguingly 
close to the value implied if collective self-gravity were the only available force to assemble 
the largest bodies (Shu 1984). 

Third, if we examine the most primitive meteorites, the ordinary and carbonaceous 
chondrites, we find no evidence for a continuous range of particulates spanning a range of 
sizes from less than a micron to, say, a meter or more (say, comparable to the size of the 
entire meteorite). Instead, we find the largest mass fraction of such objects to be composed 
of chondrules: quasi-spherical, once molten, inclusions of millimeter and smaller sizes that 
give evidence of once having been intensely and briefly heated (perhaps multiple times). 
This fact suggests that sticky agglomerative events, such as those experienced by compound 
chondrules (Rubin 2000), required special transient heating events to overcome the obstacle 
that liquids are not stable thermodynamic phases for any temperature at the low pressures 
prevalent in protoplanetary disks. In any case, from the meteoritic record, the early solar 
system failed to generate by primitive processes any compact particulates in excess of a few 
centimeters in size (the largest refractory inclusions). (See Shu et al. (2001) for a promising, 
although unconventional, mechanism for producing the refractory inclusions and chondrules 
in chondritic meteorites.) Electrostatic attraction could have played a role in building looser 
aggregates if the individual particulates acquired significant levels of charge (Marshall & 
Cuzzi 2001). This effect has been seen in zero-gravity experiments with particle densities 
well above the threshold for gravitational instability. In order to be a relevant growth 
mechanism, it must be shown that tribocharging, the balance between coUisional charging 
and ion/electron discharging, yields electrostatic attraction at much lower particle densities. 

Fourth, even if a mechanism could be found to grow chondrule-sized particulates to 
meter-sized bodies, one would have to worry about the rapid inward orbital drift associated 
with gas drag that would carry such bodies from 1 AU into the protosun on a time scale of 
only ~ 10^ yr (Weidenschilling 1977). In contrast, mm- and km-sized bodies have gas-drag 
drift times in excess of 10^ yr. Only by the direct assemblage of chondrules and related 
objects into planetesimals, avoiding intermediate steps, can one prevent a rapid loss of solid 
material from the solar nebula by gas drag. 

Such a direct-assemblage mechanism exists in the gravitational instability (GI) proposal 
put forth by Goldreich & Ward (1973). A similar theory was advanced independently by 
Safronov (1969). In the Goldreich- Ward theory, particulate settling yields a subdisk of 
solids that is thin and non-dispersive enough to make overdense regions undergo runaway 
local contraction. This occurs when Toomre's criterion for axisymmetric GI in a rotating 
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disk (the nonaxisymmetric case is similar) is satisfied: 

where Q is the angular Kcplcrian rotation rate, and c and S are the velocity dispersion 
and surface density. Throughout we use "p" and "g" subscripts to refer, respectively, to 
the particle and gas components of the disk. The Goldreich-Ward instability should not be 
confused with the mechanism of Boss (2000), who considers the formation of coreless gas 
giant planets from GI of gas disks. 

Toomre's criterion (1) is equivalent (within factors of order unity) to the "Roche" limit 
which has been derived specifically for the case of stratified fiuids(Goldreich & Lynden-Bell 
1965; Sekiya 1983). Sekiya finds GI to occur when the particulate plus gas space-density, 
p — Pg + Pp, at a distance r from a star of mass exceeds a certain critical value in the 
midplane: 

p > PR = 0.62Mjr\ (2) 

At a distance r — 1 AU from the Sun, = 4 x 10"^ g/cm^, which implies a critical space 
density of rock that is roughly seven orders of magnitude less than the internal density of 
compact rock. Thus, in the stages preceding actual planetesimal formation, we may treat 
the collection of solids as an additional ideal "gas" co-mixed with the real gas of the system. 

Operating at a radius of r = 1 AU, the self-gravitating disturbance with the most 
unstable wavelength creates ~ 5 km planetesimals in ~ 10'^ yr (Youdin & Shu 2002). The 

process occurs on a time scale longer than orbital periods (~ 1 yr in the zone for terrestrial 
planet formation) because of the need to damp spin-up and random velocities during the 
contraction to planetesimal densities. But the important point remains, that by leapfrogging 
intermediate size regimes, GI avoids the rapid inspiral of meter-sized bodies. 

Unfortunately, a powerful argument has been developed against the GI scenario, which 
has led largely to its abandonment by modern workers in the field (Weidenschilling 1995). 
A review of the difficulty is necessary before before we can justify a renewed attack on the 
basic idea. 

Even in an otherwise quiescent disk, midplane turbulence may develop to stir the par- 
ticulate layer too vigorously to allow sufficient solid settling to the midplane. Without such 
settling the criterion (2) cannot be satisfied. The problem lies in the vertical shear pos- 
sessed by disks with a highly stratified vertical distribution of solid to gas. The particulate- 
dominated sub-disk, which possesses near-Keplerian rotation, revolves somewhat faster than 
the surrounding gas disk, which has non-vanishing support against the inward pull of the 
Sun from gas pressure in addition to centrifugal effects. The magnitude of the resulting 
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velocity differential, A^;^ = tjvk — rjrQ,, is proportional to 77, which roughly equals the ratio 
of thermal to kinetic energy of the gas: 

(dP/dr) , , ,0 

where P is the gas pressure and Cg is the isothermal sound speed. 

In the popular model of the minimum solar nebula (hereafter MSN, sec §2), 77 ~ 2 x 
10~^(r/AU)^/^, and Av^ ~ 50m/s at r = 1 AU (Hayashi 1981). Turbulent eddies with a 
characteristic velocity equal to the available velocity differential, ~ Av^, would then prevent 
Gl, since the Toomre Q criterion requires that the particle random velocity be much smaller: 
Cp < 7 cm/s -C Av^, for instability (Weidenschilling 1995). 

These general arguments are supported by numerical simulations (Cuzzi, Dobrovolskis, 
& Champney 1993) which calculate the steady state properties of two-phase (gas and par- 
ticulate) turbulence in the midplane of a MSN disk. Particulates with internal densities of 
normal rock and with sizes from 10 to 60 cm (which are assumed to have grown by other 
mechanisms and are moderately coupled to gas motions) acquire space densities too low for 
GI by an order of magnitude or more. Their computational methods use a mixing length 
prescription to relate diffusivity and velocity shear through an extrapolation from labora- 
tory studies of boundary layers. While quite sophisticated, this approach docs not directly 
address the mechanism by which the existence of the vertical shear generates turbulence. 

Using a hnear stability analysis, Sekiya (1998) confirmed that GI in a turbulent dust 
layer is impossible, unless the ratio of dust to gas surface densities is significantly enhanced 
over normal cosmic values. Part of the purpose of the current paper is to understand the 
physical basis of Sekiya's conclusion. But for the present, let us merely note that there exist 
several possibilities for enhancing the solid/gas ratio in protoplanetary disks, either globally 
or by local concentration. 

First, in the X-wind model (Shu, Shang, & Lee 1996), chondrules are created from 
material at the magnetically truncated inner edge of the disk, called the X-point. Solids 
and gas arc launched from the X-point in a bipolar outflow. While the gas escapes in a 
collimated jet, solids of roughly millimeter size fall back to the disk at planetary distances, 
thus increasing the disk's solid/gas ratio. Since 1/3 of all material which passes through the 
X-point is launched in an outflow, while the remaining 2/3 accretes onto the protostar, rocky 
material in the disk could be augmented by a total amount as much as 4 x 10^° g, or 30 times 
the amount of rock in a standard MSN model. Such an enhancement factor is more than 
sufficient, as we shall see, to promote GI in the sub-disk of solids. In point of fact, because 
of efficiency considerations in the manufacture of chondrules and refractory inclusions and 
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their irradiation to produce short-lived radionuclides (see, e.g., Gounelle et al. 2001), there 
are reasons to believe that an amount of rock not much larger (but perhaps a few times 
larger) than that contained in a MSN was recycled by the X-wind to the disk. We shall find 
that what is important is the ratio of solid to gas surface densities in the disk. Once we 
accept that this ratio need not be cosmic (e.g., rock to gas = 4 x 10^'^), then there exists 
no a priori theoretical objection to a revival of the Goldreich-Ward mechanism for forming 
planetesimals. 

Second, since solids tend to settle towards the midplane, the surface layers of a disk 
should become relatively gas-rich. Thus, any mechanism that removes material from the sur- 
face of a stratified disk would increase the solid to gas ratio computed in terms of vertically 
projected column densities. Possibilities for such surface removal include (1) photoevapora- 
tion, which dominates in the loosely bound outer disk (Shu, Johnstone, & HoUenbach 1993); 

(2) layered accretion, which occurs if only the surface layers of a disk are sufficiently ionized 
to support magneto- rotational turbulence (Gammie 1996); and (3) stripping by stellar winds, 
which is probably more effective near than far from the star (Hollenbach, Yorke, & Johnstone 
2000). The amounts of solid to gas enhancements achievable by these processes are difficult 
to predict, but the timescales for the dominant processes are typically < few x 10® years, 
and thus likely to be relevant to the evolution of protoplanetary disks. 

Third, gas drag can also lead to local enhancement of particulate concentrations by a 
variety of mechanisms. (1) Isotropic turbulence with a Kolmogorov spectrum concentrates 
particles in numerical (Squires & Eaton 1991) and laboratory (Fessler, Kulick, & Eaton 
1994) experiments. Extrapolation to the high Reynolds numbers of protoplanetary disks 
imphes concentration of chondrules by factors of up to 10^ (Cuzzi et al. 2001). However, 
these estimates do not take into account the redispersal of the concentrated pockets of solids 
if the turbulent eddies are intermittent and do not maintain fixed centers. (2) Similarly, disk 
vortices could concentrate chondrules as well, but they are more effective for meter sized 
bodies (de la Fuente Marcos & Barge 2001). There also remains the issue whether vortices 
will rise spontaneously in protoplanetary disks if there are no natural stirring mechanisms. 

(3) Secular instabilities associated with gas drag might concentrate particulates, even without 
self-gravity (Goodman & Pindor 2000). Goodman and Pindor assumed that gas drag acts 
collectively on a particulate "sheet" , a valid approximation when turbulent wakes overlap. 
Unfortunately, wake overlap seems unlikely unless the particles are fairly closely packed, in 
which case GI would already be effective. In a related context. Ward (1976) has shown in 
an underappreciated study that viscous drag modifies the standard GI criterion through the 
introduction of an additional instability that can occur at values of Q 3> 1. However, this 
additional instability, being secular in nature, has a much smaller growth rate than the usual 
Goldreich-Ward mechanism. (4) In §4, we develop the simplest concentration mechanism 
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for particulates: radial migration due to gas drag. 

The fundamental assumptions required for the results of this paper are the existence, 
at some epoch of the disk's evolution, of (1) relatively quiescence in the midplane regions, 
even though the surface layers may be undergoing active accretion (Gammie 1996), and (2) 
compact solids with chondrule-like properties that are well, but not perfectly, coupled to the 
gas through mutual drag. 

Under these conditions, we argue (1) that vertical shear can only induce a level of 
midplane turbulence that has limited ability to stir solids, with the critical value of the 
surface density of solids being given roughly by Ep^c ~ V^Pgl (2) that gas drag alone 
can lead to a global radial redistribution of solids so that the local surface density of solids, 
Sp, can exceed the critical value, Epc- Therefore, whether a recycling of solids occurs 
by the X-wind mechanism or a radial redistribution of solids occurs by simple gas drag, 
we conclude that the planet-forming zones of the primitive solar system can achieve the 
requisite conditions for the formation of planetesimals on a time scale comparable to the 
typical lifetime, ~ 3 x 10^ yr, that has been inferred for the disks of T Tauri stars (Haisch, 
Lada, & Lada 2001). However, the margin for success is not large, and it could be that planet 
formation is a less universal phenomenon than it has been widely touted to be, and that it 
has a much greater diversity of outcomes (including a complete failure to form any planets) 
than suspected prior to the discovery of extrasolar planets (Mayor and Qucloz 1996, Marcy 
and Butler 1998). Indeed, the very fact that planetesimal formation may involve a threshold 
phenomenon, namely the existence of a nontrivial critical surface density, Ep c ~ i^rpg, implies 
that low-met allicity systems should be much less likely to form planets than high-metallicity 
systems, a correlation which seems already to be present in the empirical literature (Gilliland 
et al. 2000, Laughhn 2000). 

This paper is organized as follows. In § 2 we present the basic properties of our disk 
models. After reviewing the techniques of Sckiya (1998) for deriving density distributions 
of well-coupled particles in §3.1, we physically interpret, in §3.2, the midplane density 
singularities which occur in these profiles as evidence that midplane turbulence can only stir 
a finite amount of material. This allows us, in § 3.3 to calculate the solid/gas enhancements 
required for GI in various model disks. We show that aerodynamic drift can concentrate 
particles of a given size radially in the disk on cosmogonically interesting time scales in § 4.1, 
and we generalize to distributions of particle sizes in §4.2. In §4.3 we evaluate whether this 
concentration mechanism provides enough enhancement to induce GI. Closing remarks are 
made in § 5. 
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2. Disk Properties 

If any accretion and its associated turbulence occurs in the disk, we assume that they 
are confined to the top and bottom surface layers of the disk (Gammie 1996). Because 
the deeper (midplane) layers of the disk are then heated only by radiation from above or 
below, we can then model the gaseous component as being in hydrostatic equilibrium with 
a vertically isothermal distribution of temperature, while the radial distributions for the 
temperature and surface density are taken, for simplicity, to have simple power-law profiles: 

T = (4) 
Eg = 1700/,^^-^' g cm-2. (5) 

The normalization occurs relative to 1 AU; i.e., we define a dimensionless radius, w = r / AU, 
and then 280/t and 1700/^ represent the gas temperature (in K) and surface density (in 
g/cw?) at 1 AU. The power-law indices, q and p, are the remaining parameters that define 
our gas disk. In the MSN model advocated by (Hayashi 1981), fx = fg — ^, P = 3/2, and 
q — 1/2; but for completeness, we shall consider broader ranges of possible models. 

We ignore details of the vertical structure of the gas disk, because our interest lies in 
the midplane region, to which particles settle on a timescale (Goldreich & Ward 1973): 

tset ~ ^ - 10^ (^) yr. (6) 

In the above, a is the particle radius and ps is the internal density of the solid material, 
typically 3 g/cm~^ for rock and 1 g/cm~^ for ice. 

We shall initially consider rock (r) and ice (i) distributions that reflect the power-law 
distributions of the gas: 

= 7.1/rW-^gcm-^ (7) 

r if T{w)>mK 

\ 23/iro-P g cm-2 if T(ro)<170K ' 

where the total surface density of sohds is: Ep = Er + S^. At cosmic proportions, /g, /r, and 
/i would all be equal to each other (and equal to 1 in the Hayashi model), but we shall relax 
this restrictive assumption in what follows. 

Table 1 shows the parameters for the MSN and other models used in this paper. Values 
of the reference state at solar abundances are shown. The flrst letter of a model: H, A, or B; 
indicates its temperature proflle: warm (Hayashi's MSN values), cool, or cold; respectively. 
The surface densities are consistent with the MSN (p = 3/2, /r = 1) unless denoted otherwise, 
e.g. appending an "f" to indicate a flatter p—1 proflle. 
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Our temperature and surface density profiles are chosen to lie within the bounds set 
by astronomical observations, particularly the mm-wave continuum emission from T Tauri 
disks (Osterloh & Beckwith 1995), as well as with the midplane temperatures predicted by 
the theory of passive, flared disks (Chiang & Goldreich 1997). These considerations yield 
icelines considerably interior to Jupiter's 5.2 AU orbit. Uncertainties in opacities mean that 
while the particle mass in most T Tauri disks is at least that of the MSN, it could be larger 
if mass is hidden in particulates much larger than the wavelength of the observations. While 
ISO detected warm H2 around T Tauri stars (Thi et al. 2000), more detections with higher 
singal-to-noise are needed to give stringent constraints on gas mass. 



3. The Connection Between Solid/Gas Ratio and Gravitational Instability 

3.1. Dust Density Profiles 

We briefly review Sekiya's (1998) technique for deriving dust density profiles before 
interpreting the singular cusps which appear in these profiles. For small particles which 
are well-coupled to gas motions, i.e. when the particle stopping time is shorter than the 
dynamical and eddy turnover times, the gas-solid mixture can be thought of as a single 
stratified fiuid. This limit applies very well to mm-sized and smaller particles, however the 
particles must be large enough so that their settling times, see (6), are shorter than the time 
since global turbulence becomes weak enough to allow particulate settling. The well-coupled 
limit is the most conservative assumption in which to demonstrate GI since larger, decoupled 
solids are stirred less efficiently by gas turbulence. 

Radial hydrostatic balance yields an orbital velocity of the combined fluid which depends 
on the vertical distribution or particles (Nakagawa, Sekiya, & Hayashi 1986): 



VK,. (9) 



To derive the above, we have ignored the variation of the gas density in the vertical direction 
since the particulate subdisk is much thinner than a gas scale height. When solids provide 
most of the inertia, Pp ^ Pg, the fluid motion is Keplerian, but when gas dominates we have 
the usual formula for pressure-supported rotation, = (1 — r))vK- 

As particles settle towards the midplane, the vertical shear rate (which promotes mixing) 
and the buoyancy (which stabilizes against mixing) both increase. This competition between 
destabilizing and stabilizing influences is conventionally characterized by the Richardson 
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number: 

In the above formulation, the Brunt- Vaisala frequency, N , is a measure of buoyancy: 

N'' = g,d\up/dz, (11) 

where Qz is the vertical gravity. When the Ri drops below a critical value, Ric, typically 
^ 1/4, we have Kelvin-Helmholtz instability (KHI) and the onset of midplane turbulence. 

At low particulate densities pp <^ pg, we may approximate the vertical gravity as coming 
solely from the central star, ~ Vt^z, while differentiation of equation (9) yields for equation 
(10): 

Ri oc (12) 

dp^/dz 

Increasing stratification is then destabilizing, since the increase in shear outweighs the in- 
creasing buoyancy, decreasing Ri until it reaches Ric- 

In our detailed calculations, we include the self-gravity of the combined fluid in the 
plane-parallel approximation: 

= Q'^z + 47rG' / [pg + p^{z')]dz' . (13) 

When pp becomes large because of the settling of particulates to the midplane, the self-gravity 
of sub-disk can provide a significant stabilizing infiucncc against KHI. Indeed, we shall find 
that this additional effect is the cause of a sudden "precipitation phenomenon" when the 
surface density of solids crosses a critical value. Since the self-gravity of the subdisk of solids 
is also of crucial importance for the Goldreich-Ward mechanism, we find that "precipitation" 
is virtually synonymous with the satisfaction of the criterion for GI, equation (2). 

To see this, note that for the second term on the right-hand side of equation (13) to 
become of comparable importance as the first, we require 

47^(pg + pp)~QVG' = M,/r^ (14) 

which differs from equation (2) only by a factor of 7.8. Precipitation is formally reached 
before GI (if we think of a process of slowly increasing the local sohd/gas ratio in the disk), 
but since precipitation results, as we shall see, in the settling out of a very dense midplane- 
layer of solids, GI follows almost immediately after conditions in the disk become ripe for 
the beginning of precipitation. 

Before these dramatic events occur, the particulate layer should reach an equilibrium 
which is marginally unstable to the KHI, because even weak turbulence is sufficient to halt 
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the settling of well-coupled particles. We henceforth assume with Sekiya that the subdisk of 
solids satisfies at each value of the marginally stable condition, Ri = Ri^ = 1/4, consistent 
with the necessary criterion for instability in plane-parallel flows. The result is probably 
modified by the disk geometry and Coriolis forces. However, preliminary investigations of an 
analogous two-layer model (Youdin, unpublished) indicate that instability is not significantly 
altered by the inclusion of rotation, though the analog of Ric decreases slightly. In any event, 
meaningful quantities (particle scale height and critical surface density) scale as y/RQ, so 
the exact value used for the critical Richardson number does not have a strong effect on the 
numerical results. 

The ansatz of a marginally stable state, Ri = 1/4, defines via (10) and (13) an integro- 
differential equation for particle density profiles, Pp{z), at a given radius. The solution [see 
eqn. (18) in Sekiya (1998)], requires two boundary conditions: the midplane particulate 
density, Pp(0), and reflection symmetry across the midplane, dpp/dz — 0. The former is 
equivalent, as a one-to-one mapping, to an integral constraint on Ep. 

Density proflles of this type can be analyzed for GI, e.g. using (2). One typically finds 
that low mass disks are stable to the Goldreich-Ward mechanism at cosmic abundances for 
solids/gas. However increasing the surface density in particles, Ep, while holding Eg fixed, 
leads to the development of a density cusp (see Fig. 1), which quickly leads to GI in the 
midplane, as we described above. Decreasing Eg relative to Ep has the same effect. In this 
context, it is well to recall that the evidence of the planets of the solar system tells us only 
what the minimum content of the solids was in the primitive solar nebula; it says almost 
nothing about either the maximum content of solids or actual content of gas at the time of 
the formation of planets and planetesimals. 

Sekiya (1998) was the first to find the appearance of density cusps in the above type 
of analysis, but he curiously dismisses the "seeming infinite density" as "due to [an] over- 
simplified analysis." In contrast, we believe, as already hinted above, that the appearance 
of density cusps, i.e., the phenomenon of "gravitational precipitation" is intimately tied to 
the process of GI, and is indeed crucial to understanding how the formation of planetesimals 
might occur in nature. 

3.2. Saturating the Peirticulate Layer - Qualitative Understanding 

The cusps which appear in constant Ri particle density profiles become singular, i.e. 
reach infinite midplane volume-densities, for finite particulate surface-densities. Critical 
profiles cannot be constructed for higher values of Ep. We argue that the particulate layer 
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becomes "saturated" at this critical surface density and excess solids will precipitate to the 
midplane and undergo GI. 

The value of this critical surface density is found by taking the Pp(0) oo limit of 
the general result for the surface density as a function of midplane density, equation (22) in 
Sekiya (1998): 

Ep,c = 2-\/^77rpg • s(V'), (15) 

where 

s{^^J) = {l + ^^J)\n[{l + '^P + ^/T+2^)/^^J]-^/l + 2^, (16) 

is an order unity term which depends only on the (typically weak) self-gravity of the gas: 
t/j = AnGpg/^ll^ ^ 1.9 /Qy. The height of this critical layer can be found similarly from 
equation (21) of Sekiya (1998): 



Hp,c = y/W^rir-hi^), (17) 

h{ip) = ^/lT2i^-1pH{l + 1p + y/lT2i^)/iJ]. (18) 

These self gravitational terms can be approximated over a wide range of Qg values by: 
s(Qg) ^ 1.8(gg/10)°-35 and h{Qg) ^ 0M{QJ10)'^-'\ 

When the self-gravitational terms are neglected, we see that the particulate mass Sp.c ~ 
PgHp c that can be stirred by midplane shear is equal to the mass of gas in the layer. This 
result has a simple intuitive interpretation via a dust-storm analogy. A fierce wind-storm 
in the desert can pick up a lot of dust via the Kelvin-Helmholtz instabilities that ruffle the 
interface between the air and the desert floor. However, the wind cannot pick up the sand of 
the whole desert. There is a maximum amount of dust with which air can be laden before as 
much dust falls out of the air as is picked up by it. In the case of a desert storm, this amount 
depends on hard the wind is blowing. In the nebular disk, the strength of the "wind" is fixed 
by relative mechanical balance considerations (depending on the parameter rj), and it should 
not be too surprising, given the lack of intrinsic scales in the problem (when we ignore a role 
for self-gravity) , that this saturation level is roughly reached when the mass of dust is equal 
to the mass of the layer of air in which the dust is embedded. 

To derive this result more quantitatively (but still roughly), we consider how the KHI 
depends on Sp and the particle layer thickness. Hp, for low and high values of Pp/pg. We 
do this by making the crude approximation dpp/dz ~ Pp/ Hp pa Ep/i7p, which allows us to 
express: 

^ ifA+y, (19) 

When pp <^ pg, the Richardson number scales as i?i oc ifp/Ep. Thus, maintaining the 
balance between buoyancy and shear as Ep increases requires the particle layer to become 
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slightly thicker (as is seen in the detailed solutions), and will lead to a density increase: 
Pp ~ Sp/i?p oc Ep''^. 

When the particle density becomes large, Pp ^ Pg, it is no longer possible to maintain 
KHI since Ri ocT,"^. While the buoyancy is relatively constant (as long as we can still ignore 
self-gravity), the shear actually decreases with added mass, dV(f,/dz oc {dpp/dz)p~'^ oc E"-*^, 
because the velocity contrast is diminished as more of the material rotates at speeds closer 
to Keplerian. From this simplified analysis one finds that the maximum Sp which can be 
supported by the KHI is: Sp^c ~ \/RicWPg,i agreement with the detailed result except for 
the self-gravitational correction factor. 

As already hinted upon, saturation does not lead to GI of the entire particle layer 
unless the gas disk is already self-gravitating: Qp ~ Vt^H^^d (^rCSp^c) ~ Qg- Thus, only the 
unstirred particles in excess of Sp c should undergo GI initially. Self-gravity is important 
largely for the s{%l]) factor. If we were to take 'i/^ — > 0, which can be shown to be equivalent to 
setting Qz = Vt^z at the outset, then as a high density cusp develops, Pp(0) ^ pg, one finds 
that Ep oc ln(pp(0)). Thus, an infinite density cusp no longer corresponds to a finite Ep^c- 
Ignoring self-gravity entirely would have caused us to miss the saturation effect in Sekiya's 
detailed solutions. 

The fact that midplane shear can only stir a finite amount of solids is only relevant if 
the saturation point can be reached. Let us compare Ep^c to the surface density available at 
cosmic abundances : 

^r.IlP^^(^\ (20) 

where we use Eg ~ 2.4pgrCg/fK and t] ~ 1.6(cg/fK)^. Saturation requires that the thinness 
of the particulate layer, Hp/Hg ~ Cg/fx ~ 1/30, equalize the space densities (p) of particles 
and gas. Our estimate indicates that saturation may be possible at cosmic abundances (i.e. 
Ep c < Ep in the outer solar system, where the inclusion of ices yields a gas/solid mass 
ratio, (Eg/Ep)0 ~ 57. Saturation at cosmic abundances is unlikely if only rocks are present, 
(Eg/Ep)© ~ 240. To obtain more accurate answers, detailed solutions must be used to 
determine whether and how much enhancement is required. The simplified treatment shows 
the general trends that colder disks require less enhancement. It also tells us that total disk 
mass is relatively unimportant, until the self-gravity factor s{ip) comes into play. 

3.3. Required Enhancement Values 

In this section we compute more accurate values for the enhancement of solids, or the 
depletion of gas, required for the onset of GI in various disk models. For the reasons described 
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earlier, we adopt the "saturation" threshold (required to precipitate out a midplane layer of 
solids with a formally infinite space density): 

Sp > Sp,c, (21) 

as the criterion for the onset of GI. This is clearly a more conservative approach than 
considering the Roche stability criterion which only requires a high but finite midplane space 
density, equation (2). In practice, because saturation is reached so quickly once conditions 
become appropriate for precipitation, the two criteria arc virtTially identical, except for the 
cases when the gas is a considerable aid to Gl. In such extreme conditions of massive or very 
cold disks, we revert to the Toomre criterion, Qp < 1, to assess the possibility for Gl. 

A disk model which is gravitationally stable at cosmic abundances can be made unstable 
in several ways, as shown for the MSN in Fig. 2. Holding the gas content fixed, increasing 
Ep will yield Gl at a value, Sp^u, which is typically ^p,c- The amount of enhancement 
£ = Sp u/Sp needed for GI is shown for various models in Fig. 3. 

The inverse process, holding the particle content fixed and lowering Eg until GI occurs 
at a value Eg u, works equally well. The depletion factors, V = Eg_0/Eg^u, defined so that 
T> >1, are plotted in Fig. 4. 

One can think of the gas depletion or solid enhancement scenarios as refiecting the 
mechanism which gives rise to solid/gas ratio enhancements and the physical conditions 
at the time of planctcsimal formation. There is no essential difference between the two 
procedures except that enhancing solids yields a higher surface density. The reason that 
required depletion factors are larger, V > £, typically by factors of 1.5-2, is the effect of 
self- gravity represented by the s{ip) factor in Ep c. 

We find that GI requires augmenting the particle to gas ratio by factors of two to tens 
above cosmic, depending on the disk model and radial location. As expected, colder disks 
with less pressure support (and thus less vertical shear) need less enhancement. Higher mass 
disks also require smaller enhancement factors, but the effect is weaker. We find that GI 
without enhancement is possible only in the outer regions (> 10 AU) of cool disks which are 
10 to 15 times more massive than the MSN in these regions. 

Two caveats exist in the interpretation of these results. If only rocky materials, e.g. 
chondrules, are enhanced and not ices, then the fractional enhancement of solids is actually 
(1 -I- W)£ — W, where W — 3.2 (or 0) is the cosmic ice to rock ratio outside (or inside) 
the iceline. Also, if planctcsimal formation occurs in one of the high Ep scenarios, the total 
amount of solids need not exceed the MSN if the enhancement is local. 
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4. 



Drift Induced Enhancement 



Here we show that gas drag in a laminar disk causes global redistribution and concen- 
tration of small solids as they inspiral. It is surprising that this straightforward and robust 
effect has not been studied before, but the more complex case of the global evolution of 
solids in turbulent disks has been studied in numerical simulations (Stepinski & Valageas 
1996). In our case, the concentration of mm-sized sohds occurs on ~ 10^ year timescales, 
implying that chondrules could play a crucial role in triggering planetesimal formation and 
in the observed disappearance of dust disks around T Tauri stars. 

For particles with a radius a < 9A/4, where A is the gas mean free path, gas drag follows 
Epstein's law (Weidenschilling 1977). For the MSN, 9A/4 ~ .6ro^^/^ cm, so particles up to 
chondrule sizes can safely be treated with Epstein drag for all but the innermost regions. 
The drift speed due to Epstein drag in a hydrostatic gas disk is: 



where tst = Pso/(pgCg) <^ 1/Q is the stopping time. The numerical values in (22) apply to 
the MSN; more generally the radial dependence goes as v^t oc r'^ where: 



where we again specify p as the surface density powerlaw of the gas disk: Eg oc r~^. Note 
that d depends only on gas properties, which we assume to be time constant, and not on the 
evolving surface density of the solids. 

We will set the inertial factor pg/ p ^ 1 in (22) because this approximation simplifies 
the mathematics at no significant cost in realism. Since GI would occur (by the saturation 
mechanism of §3.2) if > pg, the procedure represents a factor ~ 2 error at worst in 
non-critical circumstances. We assume that any variation in drift speeds associated with 
variations in pg/p is hidden by the larger spread that occurs when we have a spectrum of 
particle sizes. With this simplifying assumption, v^t depends on the solid density and size 
of the particles, but not on their surface or space density. 



An axisymmetric distribution of uniformly sized particles with surface density S(r, t) 
(we drop the "p" subscript here) evolves according to the continuity equation: 




(22) 



d^p-q + 1/2, 



(23) 



4.1. Evolution of a Single Pcirticle Size 



(9E (9S E (9 



(24) 
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subject to an initial value for E(r, 0). This linear, first order PDE, a Cauchy problem, can 
be solved by the method of characteristics to yield a general solution: 

E{r,t)=r-''-'g{n{r,t)), (25) 

where 



n{r,t) 



1 

' d-1 



(26) 



is the initial location of a particle which winds up at radius r at time t. Note that the 
material at a fixed position r is arriving from increasingly distant locations ri with time, 
regardless of the sign oi d — 1? The function g{r) = r'^+^S(r, 0) is determined by the initial 
conditions. 

Consider a surface density which initially has a power-law profile with a cutoff at some 
outer radius: S(r, 0) = Eor"" if r < To and E(r, 0) = otherwise. The surface density 
evolves as: 

E(r,t) = Sor-'^-Vf+i-"(r,t), (27) 

if ri(r, t) < To, i.e. if the material comes from within the disk cutoff, and E(r, t) = 
otherwise. 

The concentration of disk material at a given location is: 

3/2-q+(p-n) 



ri(r,t) 



(28) 



where the requirement that ri{r,t) < still holds, or else C = 0. We have expressed d using 
(23). If the solids and gas share the same initial powerlaw distribution, then concentration 
depends only on the conditions that 3/2 > q (nearly inevitable) and that the outer edge of 
the disk has not passed, i.e. ri(r, t) < Tq. The maximum enhancement at any location takes 
the simple form: C^ax — {ro/r)^^"^"''- This can reach factors of lO's or lOO's for a disk with 
To ~ few X 100 AU, hke observed T-Tauri disks. Note, however, that the largest concentration 
factors in the inner disk might not be reached because of planetesimal formation that occurs 
as the small solids are drifting in. 

The time it takes to reach a certain enhancement level C < Cmax is: 

d-l 

\l — C 3(2-3) Ir 



^For the special case d= 1, the general solution (25) is still valid and n = r exp[vdi{r)t/r]. 
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For mm-sized and larger particles this timescale is shorter that the typical disk lifetime of a 
few 10^ years. Thus, significant concentration can be expected for chondrules and refractory 
inclusions. Smaller particles, such as the matrix material in chondritic meteorites, would not 
have enough time to be concentrated by this mechanism. 

The time evolution of uniformly sized {a — 1 mm) particles with an initial surface 
density: 



chosen to smooth the abrupt cutoff, is plotted in Figs. 5 and 6. The particles are embedded 
in untruncated gas disks corresponding to models H and Af, respectively, and n = p is 
assumed. For larger d values, which correspond to high drift speeds in the outer disk regions, 
concentration occurs on a shorter timescale. Also, when d > 1 the enhancement is most 
prominent where r[{r, t) = r^, i.e. the radius r to where the outer edge has drifted by time 
t. Thus, the profile becomes flatter and eventually inverted, resulting in ringlike structures. 
However when d <1 the inner regions are preferentially concentrated, steepening the surface 
density profiles. For the case d — 1 (not shown) the enhancement is uniform with radius (for 
r such that ri(r, t) < Tq). 

To summarize, the magnitude by which particles of a given size are concentrated depends 
on the radial extent of the disk, since this controls the amount of material available to pile 
up. The timescale depends on the drift rates, which are highest for gas that is warm and 
low density and particles that are large and compact. 



Generalizing our results to the more realistic case of a spectrum of particle sizes is 
straightforward since particles of a given size evolve independently at their respective drift 
speed if we adopt the basic premise that particle collisions do not lead to particulate growth. 
Initially we assume a distribution of particle sizes such that the surface density of particles 
with a size between a and a-\-da is aQ{r)N{a)a^da. The size distribution, initially independent 
of disk radius, is typically assumed to be a simple powerlaw, N{a) oc a^'^, with upper and 
lower size cutoffs, Omax and Omin- The mass distribution is tilted towards large (small) 
particles if s < 4 (s > 4) . 

Particles of a given size evolve according to (25), and integration over the size distribu- 
tion yields the total surface density: 




(30) 



4.2. 



Peirticle Size Distributions 




(31) 
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where we must now take into account the dependence of Ri on particle size via the drift 
speed. Care must be taken in evaluating this integral since the upper size cutoff, Omax; varies 
with R and t, because of the requirement that material not come from beyond a certain 
radius. Even for the case of a disk with no sharp edge, eqn. (30), a cutoff must be placed at 
some large finite disk radius because, when d > 1, material formally drifts in from infinity 
in a finite time. This effect is not physical since the formula for the drift speed must be 
modified when it exceeds the sound speed, which occurs beyond 500 AU for millimeter sized 
solids. An insignificant amount of material lies at such radii, and we may eliminate the 
mathematical problem by imposing a cutoff. Our final results are insensitive to the location 
of the cutoff as long as it is imposed in the exponential tail, r > Tq. 

The numerical integrations were performed using a fifth order Romberg method to deal 
with the large derivatives present in the kernel. Figs. 7 and 8 show the time evolution of 
the enhancement E(r, i)/E(r, 0) of particles with a distribution of sizes and masses that is 
characteristic of chondrules: amin = -01 mm, amax = 1 mm, s = 3. We again use (30) for 
S(r, 0), with To = 250 AU, with the gas disks following models H and Af. The main effect of 
introducing a size distribution is to reduce and broaden the amount of concentration since 
particles of different sizes drift at different rates. 

4.3. Saturation Criterion 

In order to determine whether the concentration due to gas drag only is enough to yield 
GI according to the saturation mechanism of § 3 we need to make choices about an "initial" 
size distribution of particles. If the particles are uniformly sized (a false assumption) then 
enhancement factors are much larger than required for GI. Here we make a more reasonable 
assumption that initially the disk contains solar abundances of matrix material with a size 
distribution: 0.1 /im < a < 0.1 mm with s = 4, augmented by 50% (in mass) with chondrules 
having the size distribution as above (in §4.2). The time evolution of the enhancement due 
to these two components in model Af at 1 AU is shown in Fig. 9. The maximum factor of 
10 in enhancement is almost exactly the factor 10 that is required for marginal GI inside 
the icclinc, sec Fig. 3. The resultant 5 to 1 ratio of chondrules to matrix material is slightly 
higher than the 4 to 1 ratio seen in ordinary chondrites. On the other hand, much of 
the matrix material in chondritic meteorites may be chondrule fragments, so the empirical 
enhancement may be considerably higher than represented by the nominal chondrule: matrix 
ratio of 4:1. In such a picture, planetesimal formation is triggered by the appearance (plus 
recycling and concentration) of chondrules in the nebular disk, with the ratio of chondrule 
to pristine matrix required to trigger GI larger in the inner solar system than the outer. 
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consistent with the trends running from ordinary chondrites to carbonaceous chondrutes. 

For model H at 1 AU, concentration under the assumptions described in the previous 
paragraph yields enhancement by only a factor of 7, not quite enough for GI. Note that if 
the introduction of chondrules (say, by the rock recycling mechanism of the x-wind theory) 
occurs not all at once in the beginning, but is delayed in time or is continuous in operation, 
then the enhancement factors of mm-sized bodies relative to the pristine matrix and gas 
(which drains continuously into the central star) may become significantly larger than the 
simple models illustrated above. 

At larger radii where ice is present, concentration factors depend on the size distribution 
of the icy material, a highly uncertain quantity. If the ice is not much larger than centimeter 
sized, then its lower internal density will cause it to concentrate on similar timescales as 
chondrules, and GI becomes possible. 

For the somewhat arbitrary parameters we have chosen, we conclude that aerodynamic 
drift provides a significant amount of concentration, but not definitively enough to cause 
midplane GI. Thus, in the case of our own solar system, one of the other solid/gas ratio 
enhancing mechanisms discussed in the introduction may be necessary. 

5. Discussion 

In this paper wc have shown that planetesimals can form by midplane gravitational 
instabilities despite Kclvin-Hclmholtz stirring if the ratio of particle to gas surface densities 
is increased above cosmic abundances. Wc have also presented a simple drift mechanism 
which can provide most or all of the particle enhancement required for this to occur. This 
mechanism has the advantage that it must operate in passive protoplanetary disks, and does 
not depend on assumptions about accretion physics. There are some attractive features of 
this scenario not yet discussed. 

Curiously, the mass in our planetary system appears to be significantly truncated outside 
of 40 AU, i.e. in the Kuiper belt region and beyond (Trujillo & Brown 2001), whereas T-Tauri 
disks typically extend to several hundred AU. The size disparity is even greater if Uranus and 
Neptune migrated from a location interior to Saturn's orbit (Thommes, Duncan, & Levison 
1999). The process of drift induced enhancement offers an explanation. The outer disk is 
drained relatively of its solid resources by inward particulate drift, and thereby becomes or 
remains inhospitable to the formation of planetesimals. 

Another troubhng aspect of planet formation theories is the "Type-I" migration of 
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Earth-sized and larger (but not large enough to open a gap) bodies due to density waves 
torques exerted on the gas disk (Ward 1997). When the resonant torques are assumed to 
damp locally in the disk, Earth-mass bodies (at 5 AU) migrate inwards in 10^ years in a MSN 
disk. This timescale is inversely proportional to the mass of the body, but increases sharply 
once the gap-opening mass, typically 10 — IOOM0 is reached. However the drift speed is 
proportional to Eg. Thus if planetesimals form due to depletion of gas below MSN values (or 
equivalently the enhancement of solids in a very low mass disk) , then subsequent earth-mass 
cores suffer less migration. Gas depletion by one order of magnitude would significantly 
increase the survival odds of a nascent planetary system, and still leave enough gas for the 
formation of giant planet atmospheres. 

To summarize, in conventional cosmogonies with unit sticking probabilities and cosmic 
abundances of sohds and gas in a MSN (Lissauer 1993), planetesimal formation is easy, 
occurring on a time scale of ~ 10^ yr or less, while giant planet formation is hard, requiring 
time scales in excess of the typical lifetimes of T Tauri disks, ~ 3 x 10^ yr. Overall surface 
density enhancements (gas and dust) above MSN values can speed up post-planetesimal 
growth (Thommes, Duncan, & Levison 2002), but this solution would exacerbate the problem 
of the Type-1 migration of planetary embryos (Ward 1997). It would also require a finely 
tuned mechanism to remove the considerably greater amount of extra gas and solids from 
the solar system. 

In an unconventional cosmogony, where sticking probabilities arc zero (except for the 
special mechanisms that produce chondrulcs and refractory inclusions near the protosun), 
where gas is depleted, and solids are enhanced relative to standard MSN values, the condi- 
tions for the formation of planetesimals and giant planets might be intimately tied to the 
gas-dust evolution of the nebular disks of T Tauri stars. The two processes would then 
naturally acquire similar time scales, related by a single continuous process of gravitational 
growth in a gas-dust disk. Attractive byproducts of this unconventional approach would be 
a corresponding alleviation of the problems of Type I migration and gas-disk dispersal, as 
well as a possible understanding of why the Kuiper belt marks a sudden apparent truncation 
of the primitive solar system. 
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pendix and other improvements. Financial support was provided by a grant from the NASA 
Origins of Solar Systems Program. A. Y. acknowledges support from an NSF Graduate 
Fellowship. 



APPENDIX 



A. CoUisional Effects 



In most of this paper we have ignored the effects of colhsions. Here we argue that this 
approach is is dynamically justified. The time between collisions is: 

ten - - ^. (Al) 

which agrees within factors of order unity with the more exact results presented in Wetherill 
& Stewart (1993) if we use Cp/ Hp Q in the second equality Thus the coUisional timescale 
is longer than the stopping time, icoii/4t ~ ^g/^p by at least an order of magnitude. 

Even if coUisional dynamics were able to introduce an effective viscosity, i/cou — Cptc 
(Goldreich & Tremaine 1978), the characteristic diffusion velocity, Vcou ~ Vcowl"^ is less than 
the drift velocity due to Epstein drag: 

— ^ r^l^ < 10-^. (A2) 

Vdr ^p 

If collisions could change the size distribution of solids due to fragmentation or mergers 
then our results for drift induced enhancement (§ 4) would have to be modified by including 
a coagulation equation, e.g. Kenyon & Bromley (2002). However we already argued in the 
introduction against the effectiveness of coUisional agglomeration. As for fragmentation, 
the presence of significant numbers of intact chondrules (the species most important for the 
enhancement) in meteorites shows that many were able to survive their coUisional history in 
the nebular disk without much shattering. 
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Fig. 1. — Vertical profiles of particle density at 1 AU, for small solids stirred by the Kelvin- 
Helmholtz instability. The gas is described by model A and held fixed. As Sp increases 
from 30 g/cm^, the value at solar abundances (with ice), a midplane density cusp develops, 
becoming infinite at Ep = Ep^c = 94 g/cw? for this model. Also shown are the gas, pg, and 
the Roche, pr, densities. 
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Fig. 2. — The MSN disk model can be made gravitationally unstable by either enhancing the 
solids (increasing Ep as indicated by the dashed arrows), by depleting the gas (decreasing 
Sg as indicated by the solid arrows), or by some combination of the two (not shown). Note 
that instability is local, and need not occur simultaneously at different radii. 
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Fig. 3. — Enhancement factor, S, required for GI vs. radius for various models. The 
discontinuity occurs because larger enhancements are needed inside the iceline. Kinks in 
higher mass models (Af,Af5) occur when Qp < 1 occurs before saturation. 
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Fig. 5. — Time evolution of the surface density of mm-sized particles with sohd densities, 

Ps = 3 g/cm^, due to Epstein drag. The gas is assumed fixed at MSN values. Since 
d = 3/2 > 1, the Sp profiles, initially described by eq.(30) with — 200 AU, become 
flatter and eventually inverted into ringlike structures. 
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Fig. 6. — Same as Fig. 5 except the gas is described by model Af. Since d = .87 < 1, Ep 
evolves to steeper profiles. 




Fig. 7. — Time evolution of the surface density enhancement of soUd particles with a size 
distribution characteristic of chondrules. The gas disk is described by the MSN. See text for 
details. 
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Fig. 8. — Same as Fig. 7 except the gas is described by model Af. 
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Fig. 9. — Enhancement of solids vs. time at 1 AU in model Af for a particle disk that is 1/3 
chondrules and 2/3 matrix material by mass. 
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Table 1. Model Parameters 



Model 


P 


r 


<1 


Ti (K) 


iceline (AU) 




H 


3/2 


1 


1/2 


280 


2.7 


1.5 


A 


3/2 


1 


.63 


170 


1 


1.37 


B 


3/2 


1 


3/4 


100 


.5 


1.25 


Af 


1 


1 


.63 


170 


1 


.87 


Af5 


1 


5 


.63 


170 


1 


.87 



— fg — fi — fi for solar abundances. 
^See §4, eq. (23). 



